Magnetization Process of the Classical Heisenberg Model 
on the Shastry-Sutherland Lattice 
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We investigate classical Heisenberg spins on the Shastry-Sutherland lattice and under an external 
magnetic field. A detailed study is carried out both analytically and numerically by means of 
classical Monte-Carlo simulations. Magnetization pseudo-plateaux are observed around 1/3 of the 
saturation magnetization for a range of values of the magnetic couplings. We show that the existence 
of the pseudo-plateau is due to an entropic selection of a particular collinear state. A phase diagram 
that shows the domains of existence of those pseudo-plateaux in the {h, T) plane is obtained. 

PACS numbers: 75.10.Hk, 75.30.Kz, 75.40.Cx, 75.60.Ej 



I. INTRODUCTION 

A spin system is frustrated when all local interactions 
between spin pairs cannot be satisfied at the same time. 
Frustration can arise from competing interactions or/and 
from a particular geometry of the lattice, as seen in the 
triangular lattice. Frustrated systems have been inten- 
sively studied over the last decades and they were found 
to present very rich behavior such as large ground-state 
degeneracies. A recent review on frustrated systems is 
presented in Ref . [l|. 

The Shastry-Sutherland lattice (SSL) was considered 
more than 20 years ago by Shastry and Sutherland as 
an interesting example of a frustrated quantum spin sys- 
tem with an exact ground statci^ It can be described as 
a square lattice with J' antiferromagnctic couplings be- 
tween nearest neighbors and additional J antiferromag- 
nctic couplings between next-nearest neighbors in every 
second square (see Fig. [1] left panel). 
This lattice attracted much attention after an experimen- 
tal realization was identified in SrCu2(B03)2»^ a com- 
pound which was first synthesized by Smith and Keszlerii 
SrCu2 (603)2 exhibits a layered structure of Cu(B03) 
planes separated by magnetically inert Sr atoms. The 
Cu^"*" ions carry a spin 5 = 1/2 and are located on a 
lattice which is topologically equivalent to the SSL (see 
Fig. [1] right panel). The quantum SSL was widely stud- 
ied both theoretically and experimentally (see Ref. for 
a review). As the ratio J' / J varies, it presents a rich 
phase diagram with quantum phase transitions between 
the gapped dimcr singlet ground state (J'/ J < 0.68) that 
was originally discussed by Shastry and SuthcrlandfS, a 
plaquette resonating valence bond state {J' / J < 0.86), 
and a gapless magnetic statd^ii Theoretical studies on 
the quantum SSL predicted the existence of plateaux in 
the magnetization curve for the rational values of the 
saturated magnetization M/Mgut = 1/8, 1/4, 1/3 and 
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FIG. 1: (color online) The Shastry-Sutherland lattice (left 
panel) and the topologically identical structure realized in 
the (001) plane of SrCu2 (603)2 and rare-earth tetraborides 
(right panel). J' bonds are the magnetic couplings along the 
edges of the squares (clear dashed lines) and J bonds (black 
bold lines) are the diagonal dimer couplings. 



1/2: Momoi and TotsukaSiS showed that the stripe or- 
der of the dimer triplets explains the M/Mgat ~ 1/2 
and 1/3 plateaux and the latter was predicted to be the 
broadest. Experimentally, magnetization plateaux were 
observed^i^ for M/M^at = 1/8, 1/4 and 1/3 (see, how- 
ever, Refs. [l^fisl for recent experimental controversies). 
Miyahara and Uedai^ explained the existence of all those 
plateaux as a consequence of the crystallization of excited 
triplets, although e.g. the existence of the M/Msat = 1/8 
plateau is not confirmed in a recent investigation^^ 

Recently, rare-earth tetraborides RB4 have been at- 
tracting new interest to the SSL,i^ii2ii^iiSi2ai^i2^ These 
compounds present a large total magnetic moment > 
1, which justifies a classical model although magnetic 
anisotropics are likely to be relevant for a description of 
the experiments. RB4 crystallize in the tetragonal struc- 
ture with space group PA/mbm. R ions are placed in 
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the (001) plane on a sublattice which consists of R-R 
dimers that are alternatively orthogonal along the [110] 
axis. This sublattice is again topologically equivalent to 
the SSL (see Fig. [1] right panel). 

Magnetization plateaux were found for various ratio- 
nal values of the saturated magnetization. Due to its 
large total magnetic moment {J = 6) TbB4 can be 
considered as a classical system. Several magnetization 
plateaux were found, including the magnetization values 
M/M,at = 1/3 and 1/2 Ref. [13 shows that ErB4 
(with total magnetic moment J = 15/2) exhibits a mag- 
netization plateau at Af/M^at = 1/2 for both h//[001] 
and h//[100]. The compound TmB4 also presents sev- 
eral magnetization plateaux for h//[001] .^^'"'"^'^°i^"'" 
In this paper we investigate classical Heisenberg spins on 
the SSL under a magnetic field. Note that we are not 
aware of any previous study of the classical limit in the 
presence of a magnetic field, although some of the first 
studies of the zero-field properties of the SSL did start 
from this limit i^i^^ In Section [ll] we discuss the classi- 
cal ground state at zero temperature and then study the 
role of thermal fluctuations as a scenario for the appear- 
ance of a M/Msat = 1/3 plateau. Note that this plateau 
arises only for non-zero temperature and is not related to 
an energy-driven mechanism. To emphasize this differ- 
ence and the fact that this plateau is not completly flat, 
we refer to it as a pseudo-plateau in the following. Then 
in Section Hill we focus on the special ratio J'/ J =1/2 
and study by means of classical Monte-Carlo simulations 
the influence of the system size and of the temperature 
on those pseudo-plateaux. A phase diagram is obtained 
and we draw a parallel between the SSL and the classical 
triangular lattice. Finally, in Section HVl we show that a 
small modification e of the ratio of the magnetic couplings 
J' /J = 1/2 + e slightly changes the phase diagram but 
still allows a pseudo-plateau to exist at M/Mgat = 1/3. 



II. THE MODEL 

We study the SSL in the classical limit. Heisenberg 
spins are represented by simple vectors Si with ||Si|| = 1. 
The magnetic field is applied to the z— axis, h = he^ 
(i.e. h//[001] in the experiments carried out on RB4). 
The sums are taken over nearest neighbor pairs according 
to the edges of the squares with magnetic coupling J' 
and over the diagonals with diagonal couplings J. The 
Hamiltonian is given by: 



J' 
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edges 



N J N 
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(1) 



By minimizing the classical energy in the absence of ap- 
plied magnetic field the lowest energy configuration is 
found to be coplanari^ It is Neel ordered for J' /J > 1, 
else it is a spiral state with an angle ip — n i: arccos (^) 
between nearest-neighbor spins. Those two possible opti- 
mum values for ip give a discrete chiral degeneracy to each 




FIG. 2: (color online) Lowest energy configuration of the SSL 
in the (001) plane for the particular ratio J'/ J = 1/2. The 
12 spins unit cell is represented in bold clear lines. Left 
panel: The 120 ° structure with three kinds of spins orien- 
tation. Right panel: The UUD state at M/Msat = 1/3. Each 
triangle contains two spins up and one spin down. Up spins 
are represented in red (dark gray) and down spins in yellow 
(hght gray). 



triangle. The choice of the angles cp in two neighboring 
triangles determines the direction of the helix. Four he- 
lices are possible which creates a supplementary four-fold 
discrete degeneracy in addition to the continuous onei^ 



A. Classical ground state at T = 

A particular ratio of the magnetic couplings allows 
the Hamiltonian to be expressed as a sum of elemen- 
tary plaquettes j^^i^^ In the SSL these are triangles which 
share edges along the diagonal dimer bonds and other- 
wise corners. When the diagonal couplings are twice as 
big as the edge couplings (i.e. J' /J — 1/2) the Hamilto- 
nian Eq. IT]) can be written as a sum over triangles up to 
a constant term: 



(J7J=l/2), (2) 



where Sa = X^^gA t'^^ total spin of one triangle and 
A^A is the number of triangles (TVa = A^). The classical 
ground state at zero temperature and in the absence of a 
magnetic field is a coplanar configuration with an angle 
Lp = ±27r/3 between two nearest-neighbor spins. There 
are three possible spin orientations S^, Sb, and Sc and 
the unit cell contains 12 sublatticcs (see Fig. [2] left panel). 
Minimizing the energy on a single triangle, one obtains 
the constraint: 



Sa = ^ , (J7J=l/2). 



(3) 



The classical ground state is obtained when this con- 
straint is satisfied in every triangle. The saturation field 
/igat is determined by the condition 5^ = 3 which gives 
/isat = 9 J'. At this value all the spins are aligned with 



3 



*z 




Az 




FIG. 3: (color online) From left to right: the "umbrella" and 
the "Up-Up-Down" states are possible spin configurations at 
M/Mgat = 1/3. The "K— state" and the canted state are 
the quasi-long-range-ordered configurations discussed in Sec. 
ElBl 



the 2— axis ("C/i7/7 state"). We focus on the field range 
< h < hsat- According to the classical constraint Eq. 
Q the classical ground state has only the typical global 
rotation as a degeneracy. 

Let us have a closer look at the particular value 
M/Msat = 1/3. The magnetic field to apply is: 



/11/3 = 3J' , (J' /J = 1/2) 



(4) 



Stable configurations at h = hi/^ must verify Sa = e^. 
Very different spin configurations satisfy this require- 
ment: the "umbrella configuration" and the "Up-Up- 
Down" state (see Fig. [3]) . In the umbrella configuration 
the three kinds of spins raise as the field increases and 
they always have the same projection on the z— axis (at 
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Sq = 1/3). In this case, 



the picture of the classical ground state with three kinds 
of spin orientations remains imtil saturation (see Fig. [2] 
left panel). On the other hand the "Up-Up-Down" state 
[UUD state) is a collinear state in which each triangle 
contains two spins Up and one spin Down (see right panel 
of Fig. . This classical ground state is required in order 
to have a classical plateau at M/Mgat ~ 1/33. One can 
easily show that at zero temperature and M/Mgat = 1/3, 
both configurations, umbrella and UUD, have the same 
classical energy. Since energetic considerations do not fa- 
vor the UUD state its existence is restricted to the field 
value /11/3. Hence no magnetization plateau can appear 
at zero temperature. 

Minimizing the energy on triangles naturally leads to 
a comparison with the classical triangular lattice, with 
magnetic coupling Ja between nearest neighbors. This 
lattice was the subject of many theoretical studies i^iSSii^ 
Kawamura and Miyashita^^ studied classical Heisenberg 
spins on the triangular lattice in the presence of a mag- 
netic field along the z— axis. They obtained classical con- 
straints that are strictly equivalent to Eq. ^ for the SSL. 
Minimizing the energy of a single triangle they showed 
that the classical ground state is completely specified by 
the conditions: 



\\Sa\\ - llSsi = IIScll 
Sa + Ss -t- Sc 



These constraints give information on the magnetization 
process of the classical triangular lattice and SSL at zero 
temperature. The application of the magnetic field does 
not break the U{1) symmetry of the 120° configuration 
which means that the continuous degeneracy of the clas- 
sical ground state remains. In other terms, at zero tem- 
perature and under a magnetic field, the 120° structure 
raises in the umbrella configuration that closes as the field 
increases and reaches the UUU configuration at satura- 
tion. 

For non-zero temperature the scenario of the magnetiza- 
tion process is completely different. Kawamura3^ showed 
that the non-trivial degeneracy no longer persists on the 
triangular lattice. Even in the low-temperature limit, the 
most favorable spin configuration is determined not only 
by its energy but also by the density of states just above 
the ground state. This entropic selection is responsible 
for the appearance of new phases which will be detailed 
for the SSL in Section HHl 



B. Effect of thermal fluctuations 

The Monte Carlo (MC) simulations reported in Sec. lIIII 
below indicate that, at finite temperature, the UUD state 
is the favored configuration for a magnetic field range be- 
low and until h — /11/3. This is possible if thermal fiuc- 
tuations raise the entropy of the UUD state relatively to 
adjacent states. Kawamura and Miyashit a^^'^^ showed 
that in the case of the triangular lattice non-zero temper- 
ature lifts the degeneracy of the classical ground state in 
the presence of a magnetic field. Also on the SSL only a 
discrete degeneracy due to the possible chiralities of the 
triangles remains. 

The selection by thermal fluctuations of particular 
states among a degenerate manifold is called Order by 
disorder^^'^^^ This selection mechanism is known to be 
at work in various frustrated systems such as the frus- 
trated square lattice^^ the kagome lattice (^^i^i^^i^ or 
the checkerboard latticci^ 

We will see that the SSL differs in some respects from the 
aforementioned cases. It is nevertheless useful to recall 
some relevant features of the order-by-disordcr selection 
before presenting the technical details of our analytical 
study of the thermal selection of the UUD state. Under 
a small but non-zero temperature the classical spin-wave 
spectrum of the selected state has the particularity to 
present soft modes. The energy of those modes is quartic 
in spin deviations instead of quadratic. Therefore they 
contribute less than "normal" modes to the free energy 
and it is clear that a state with more soft modes will be 
favored. More precisely, let us consider a state with clas- 
sical energy Eq that has A'4 soft modes and N2 = N — N4 
"normal" modes. The free energy is given by: 



3Ja 



{h < hgat) ■ (5) 



(6) 
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Ref. |3J details how soft modes should be enumerated and 
how they affect the specific heat. In a system with no 
soft modes each degree of freedom brings a contribution 
of ks /2 to the specific heat. In the quadratic approxima- 
tion the contribution of the soft modes is zero. Taking 
into account higher orders by adding quartic contribu- 
tions each soft mode adds A:b/4. 

Wc calculated the spectrum of the UUD state in the 
quadratic approximation in thermal fluctuations by ap- 
plying small deviations on the spin coordinates from the 
collinear UUD statei^ At zero temperature the 12 sub- 
lattices of the unit cell presented in Fig. [5] (right panel) 
are collinear with the z— axis. Under fluctuations their 
new coordinates are expressed in their own frame as: 



S,(r,)-(6f(r,),enr.),l-«,(r,)) 



(7) 



where = 1/2 ((ef)^ -I- (ef)^) is verifying the condition 
II Si II = 1 up to quadratic order. 

Following Ref. [SJ the Hamiltonian is expanded in spin 
deviations from the UUD state: 



1-L = E, 



UUD 



n>2 



(8) 
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FIG. 4: (color online) Lowest branch of the fluctuation ma- 
trix ^A above the UUD state in Fourier space [k^, ky). In 
the darkest regions, the black lines represent the lines of soft 
modes whose analytical expressions are given in the text. 



where EuuD is the classical energy and 7i„ ^ 0(e"). Up 
to second order in those fluctuations the Hamiltonian 
becomes, in Fourier space: 



^2 = 



k,i?— a;,y 



V*(-k)XV,(k) 



(9) 



Fluctuations are applied to the 12 sublattices of the unit 
cell and the vectors V,*(— k) [v = x,y) read: 



V*(-k)==(6l(-k),...,6- (-k)) 



(10) 



is a 12 X 12 matrix. Fluctuations act exactly the 
same way on the x- and y-components and as a conse- 
quence the matrices are the same for fluctuations in the 
X- and y-directions. Eigenvalues should be determined 
by solving an order 12 polynomial which cannot be done 
analytically. However, we can determine the soft modes 
from the zeroes of the determinant which reads: 

det(Al) (X (—2 + icoskx — 3cos2fca; -I- cosSfc^- -I- cosky)^ 

(11) 

This result shows that we have lines of soft modes, but 
no full branch. Yet this result is still exceptional since 
usually one finds only a finite number of soft modes. 
Wc further performed numerical diagonalization of M at 
each point in the Fourier space (kx,ky). For each point 
the smallest eigenvalue was selected in order to obtain 
a picture of the ground state. Fig. 2] shows the picture 
obtained over a couple of Brillouin zones. The darkest 
regions correspond to the lowest points and we indeed 
observe lines of soft modes corresponding to a cancella- 
tion of the quadratic energy. Note that the corresponding 
zero eigenvalues are non-degenerate. The height of the 
barriers between those lines is AE « 0.1 meaning that 



low-temperature behavior is expected to appear only for 
T < 0.1. The analytical expression of the lines of soft 
modes can be obtained by finding relations between k^ 
and ky so that the determinant Eq. (jlip vanishes: 



fca; = ± arccos 



1 {l±iV3){l-u{ky)) ^ 

2 4:u{ky) 



(12) 



where u{ky) = ( -\/ cos^ fcy — 1 — cosky)^^/^ . 
Following the calculations performed by Champion and 
Holdsworth for the pyrochlorc lattice^^ the specific heat 
Ch per spin is predicted to be given by: 



Ch 



NksT 



1 



7 



(13) 



where 7 is a constant related with the number of lines of 
soft modes. 

The manifold of soft modes for the SSL with magnetic 
couplings verifying J' / J = 1/2 scales like L in a 
Fourier space. Eq. clearly shows that the drop in 
the specific heat, or in other terms the effect of the soft 
modes, will disappear in the thermodynamic limit. 
The lines of soft modes are reminiscent of the q = 
UUUD state in the frustrated square latticci^ However, 
the precise mechanism is somewhat different in the SSL 
since here the soft modes do not correspond to continu- 
ous deformations of the ground state. Nevertheless, the 
UUD state is still selected for entropic reasons, like in 
the triangular lattice i^^i^^ the soft modes simply result 
in a free energy which is lower for the UUD state than 
for non-coUinear phases. The triangular lattice with clas- 
sical Hciscnberg spins docs not have full branches of soft 
modes cither, but the manifold of soft modes consists of 
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FIG. 5: (color online) Susceptibility and magnetization (inset) 
of the SSL at J'/ J ~ 1/2 and T — 0.02 for various system 
sizes. 
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FIG. 6: (color online) Susceptibility and magnetization (inset) 
of the SSL at J' /J =1/2 and for a system size L = 18. The 
temperature T varies. 



points in Fourier space. Still, the triangular lattice ex- 
hibits pseudo-plateaux at finite temperature. 



III. MONTE CARLO SIMULATIONS 

We performed MC simulations using a standard single- 
spin flip Metropolis algorithmi^ As a small refinement, 
we propose changes of the spin projection along and per- 
pendicular to the field direction independently. Pseu- 
dorandom numbers were generated by the Mersenne- 
Twister random number generator We studied square 
samples {L^ = Ly = L = v^) with periodic boundary 
conditions. Their sizes are chosen so as to be commensu- 
rate with the 12 sublattice unit cell represented in Fig. [51 
L = 6, 12, 18, 24, and 30. In the simulations the diagonal 
coupling J is set to 1. MC simulations computed in par- 
ticular the magnetization, the susceptibility x = dM/dh 
and the specific heat Ch- AH physical quantities are nor- 
malized per spin. Data was collected with up to ^ 10^ 
MC sweeps per point depending on the size and temper- 
ature considered. In order to insure thermal equilibrium, 
a comparable amount of sweeps was discarded before col- 
lecting data. Both field and temperature scans were run. 
Error bars are determined from independent MC runs 
(the average number of simulations per point is five). 



A. Magnetization pseudo-plateau at M/Msat = 1/3 

Fig. [5] shows the susceptibility and the magnetiza- 
tion as a function of the system size at J'/ J = 1/2 for 
T = 0.02. A pseudo-plateau is observed in the magneti- 
zation curve at M/M^at = 1/3. Following the previous 
analytical discussion this means that the UUD config- 
uration is thermally selected at this temperature. Its 
existence is no longer strictly limited to one single point 
at h = /ii/3, but it exists for a range of magnetic field 



values below and until /11/3. The width of those pseudo- 
plateaux increases with system size. One clearly observes 
three peaks in the susceptibility curves which become 
sharper and higher as the system size increases. The two 
first ones indicate the presence of phase transitions to the 
UU D state while the last one corresponds to the transi- 
tion to saturation. 

The stability of the susceptibility peak positions and 
height for L > 18 show that L = 18 captures reason- 
ably well the behavior of the plateau region in the ther- 
modynamic limit. Fig. [5] shows the effect of the tem- 
perature on the susceptibility and magnetization for a 
L = 18 system. The magnetization curve should tend to 
a straight line as the temperature decreases. This is con- 
firmed by the MC simulations: at very low temperature 
(see T = 0.01 in Fig. [5]) the pseudo-plateau is extremely 
small. The susceptibility presents two very sharp peaks 
around /11/3. As the temperature increases the pseudo- 
plateau becomes broader, which confirms the entropic 
selection of the UUD state in the plateau phase. Finally 
higher temperature (see T = 0.1 in Fig. [6]) destroys it. 
The peaks in the susceptibility curve become rounded 
and then completely disappear. 

The specific heat Ch (Fig. [7]) exhibits a first peak at- 
tributed to the transition from the coUinear UUD state 
to a low magnetic field state that will be detailed in Sec- 
tion nil Bl MC simulations on system sizes L = 6 — 42 
showed that the height of this peak scales with increas- 
ing system size as ^ L""/^ with a/t/ <C 1 which suggests 
a continuous phase transition. Note that the MC simula- 
tions presented in Fig. [7] were performed at T = 0.02 
which is a too high temperature to have a convergence of 
the specific heat towards 1 as the system size increases, 
as one would expect from the equipartition theorem and 
Eq. 1131 The second peak corresponds to the transition 
to saturation. 
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FIG. 7: (color online) Specific heat of the SSL with J' /J = 
1/2 at T = 0.02 and for various system sizes. 



B. Phase diagram for the ratio J'/J = 0.5 

Next we discuss the domain of existence of the UUD 
collinear phase in the {h, T) plane. Up to 5 MC simula- 
tions per point were performed on a L = 12 (TV = 144) 
system which gives a good qualitative picture of the po- 
sition and nature of the phases. Both temperature and 
magnetic field scans were performed. The boundaries of 
the phases are determined by the positions of the peaks 
in the susceptibility x- Fig. [8] shows this schematic phase 
diagram in the plane {h,T). Numbers are attributed to 
the three susceptibility peaks observed from low to high 
magnetic field in magnetic field scans of Fig. [6] (curves 
with circles) and to the two peaks from low to high tem- 
perature in temperature scans (curves with triangles). 
We simulated by means of MC the spin texture for 
L — 12,18 systems in order to obtain a qualitative picture 
of the spin configuration in the phases below and above 
the UUD phase. Fig.[9]shows some snapshots of the spin 
configurations in a L = 12 system at T ~ 0.003. Like in 
the classical triangular lattice the magnetization process 
at finite temperature is very different from the one at 
zero temperature. The phase in the low-field region is 
the "Y configuration" (see Fig. [3]). In this configuration 
each triangle contains two spins above the (a;, y) plane 
(red dark circles in Fig. [9]) while the last one is point- 
ing down and is almost collinear with the z— axis (yellow 
light circles in Fig. [9]) . This state is characterized by a 
single angle between the two spins pointing in the positive 
z— direction. The Y configuration breaks the rotational 
symmetry around the z— axis and as a consequence of the 
Mermin- Wagner theorem^ cannot be long-range ordered 
for T > 0. The spins above the (x, y) plane arc at most 
quasi-long-range-ordered with correlation functions that 
decay algebraically. On the other hand one can expect 
long-range order on the remaining down spin sublattices. 
Then in a magnetic field range below and until hi/^ the 
collinear UUD state is the lowest energy configuration 
(see center panel of Fig. As the temperature goes to 



5t . , . , . , 

-e- 1^' peak ofx 

■sat" e — ^^ 2"''peakof j: 

4 ^ a^** peak of 

^^^N,^ A l^^peakof;t; 




0.0001 0.001 0.01 0.1 



FIG. 8; Phase diagram of the classical SSL for the ratio of 
the magnetic couplings J'/J — 1/2. 



zero, the width of this phase converges to a single point 
at exactly h = hi/^ = 3/2 which is in complete agree- 
ment with the analytical prediction. This phase does not 
break the U{\) symmetry and hence true long-range or- 
der is realized. 

In the high-field region the ground state is the "canted 
state" (see Fig. [31). In this state the spins down are raising 
as the magnetic field increases while the spins up are no 
longer collinear with the z— axis. This state is character- 
ized by two angles: one between the z— axis and the two 
spins above the {x, y) plane and one between the z— axis 
and the single spin below the (x, y) plane (see bottom 
panel of Fig. [9]). As in the Y phase, at most quasi- long- 
range order is expected at finite temperature. 

In the higher field region the system reaches satura- 
tion and all spins are pointing up coUinearly with the 
z— axis [UUU state). The zero-temperature limit of the 
saturation field is /isat = 9/2, as derived in Section fll Al 
At higher temperature the system enters the disordered 
paramagnetic phase. The highest temperature where the 
UUD state still exists T* is estimated as T* w 0.17. 
Note that the classical spin-wave spectrum contains en- 
ergy barriers between lines of soft modes of the same 
order of magnitude. Hence this temperature agrees with 
the range of temperatures in which we expected the low- 
temperature behavior to appear. 

The phase diagram of the SSL with the particular cou- 
pling ratio J'/J = 1/2 presents similarities with the 
one of the classical triangular— and kagome^ lattices. 
Watarai et al^ suggested that in the classical triangu- 
lar lattice with Hcisenbcrg spins these transitions could 
be of the second order. The scaling of the specific heat 
shown in Fig. [7] is in agreement with a continuous transi- 
tion. We propose that a Berezinskii-Kostcrlitz-Thouless 
transition^i^i^ takes place as the system enters the 
UUD phase. Following what was obtained in the case 
the triangular lattice, the transition from the collinear 
UUD state to the disordered phase should belong to the 
universality class of the three-state Potts model4^i^ 
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FIG. 9: (color online) Snapshots of the MC simulations of the 
spin texture in a L = 12 system for J'/ J = 1/2 and T ~ 0.003. 
The circles represent the projection on the axis, red (dark 
gray) if S'' > 0, yellow (light gray) otherwise. The radius is 
proportional to \S~\. The arrow denote the projection in the 
lattice plane. From upper to lower panel: y-configuration 
at h — 1.2, UUD configuration at hi/3 = 1.5 — /isat/3 and 
canted state a.t h — 1.8. 
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FIG. 10: (color online) MC results for the susceptibility and 
magnetization (inset) of the SSL at J'/ J — 0.4 and T = 0.02 
for various system sizes. 



IV. STUDY OF THE RATIO J'/ J = 1/2 + e 

In this section we consider now a small deviation e from 
the ratio J' / J = 1/2. In the following discussion we use 
Monte Carlo results for J'/ J = 0.4 (i.e. e = -0.1). We 
observed that qualitatively similar behavior also appears 
for the other sign of e. 

We performed MC simulations on systems sizes L = 
6,12,18,24, and 30. Data were collected with up to 
10^ MC steps per point. Fig. [TO] shows the susceptibil- 
ity and the magnetization as a function of the magnetic 
field at the same temperature as in Fig. [5] for the ratio 
J' /J = 1/2. The susceptibility still presents two peaks 
around M/Mgat = 1/3 and a pseudo-plateau is observed 
in the magnetization curve. Therefore one can expect 
UUD to be the selected configuration at T = 0.02 even 
for ratios J' /J = 1/2 + e. According to the position of 
the peaks of the susceptibility, the width of the pseudo- 
plateau increases as the system size increases. 
Following exactly the same calculations as the ones pre- 
sented in Section FlIBI for J' /J = 1/2, we applied ther- 
mal fluctuations on top of the UUD state for the ratio 
J' /J = 0.4. The fluctuations matrix A4 exhibits negative 
eigenvalues which means that the UUD configuration is 
no longer the favored configuration in the quadratic ap- 
proximation for the spin deviations. In other words, the 
collinear UUD state is no longer selected at "very low 
temperatures" for J'/ J 7^ 1/2. We used MC simulations 
in order to obtain a picture of the spin configuration in 
the temperature range T < 0.001. It turns out that the 
favored configuration is an umbrella that is closing un- 
til saturation as the magnetic field increases. At "higher 
temperature" (T > 0.01) we observed that for h « /11/3 
the system enters the UUD phase. This "high temper- 
ature" regime is not analytically accessible with lowest- 
order thermal fiuctuations. 

Following the same procedure as in Sec. IIII Bl wc per- 
formed magnetic field and temperature scans in order to 
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FIG. 11: (color online) Phase diagram for the SSL for the 
ratio J'/J = 0.4. 

determine the region in which the UUD state survives. 
Fig. [TT] shows a schematic phase diagram for the SSL at 
J'/J = 1/2 + e = 0.4 in the {h,T) plane. Five separate 
simulations per point were performed on L = 12 sys- 
tems. MC data were collected up to 10'' MC steps per 
point. Temperature (circles) and field (triangles) scans 
were performed with the same ordering of the suscepti- 
bility peaks as in Section [III Bl As one expects from the 
analytical arguments, the resulting phase diagram differs 
from the case J'/J = 1/2. In the lowest temperature 
region, the system undergoes the same kind of magneti- 
zation process as at zero temperature: the coplanar spiral 
becomes an umbrella configuration as the magnetic field 
increases. The angle between two nearest-neighbor spins 
is TT ± arccos (1/2 -I- e). Hence, the spiral at zero temper- 
ature is not exactly commensurate with a 12 spin unit 
cell. However, as the temperature increases the spin po- 
sitions fluctuate and a 12 sublattice unit cell is recovered 
in average on the spin positions. Hence we suggest that 
the transition from the umbrella to the phases with 12 
sublattices (squares in Fig. Ilip is an incommensurate- 
commensurate phase transition. 

At higher temperature wc recover the three phases that 
were obtained for J'/J = 1/2 in different magnetic field 
regions. In the low-field region we still obtain the "Y 
configuration" and in the high-field region, before satu- 
ration, the "canted state" . The long-range-ordered UUD 
coUinear state still survives in a small region between 
the two aforementioned phases. For T < 0.01 the UUD 
phase is expected to become narrower as the temperature 
decreases (dashed bold lines). A precise determination 
how the four phases merge (black circle) is beyond the 
scope of this work. 



V. CONCLUSION 

We studied the magnetization process of the Shastry- 
Sutherland lattice in the classical limit. We found a 
pseudo-plateau at 1/3 of the saturated magnetization 
which corresponds to a coUincar UUD state. The spec- 
trum of spin waves above this state has lines of soft 
modes, like the q = UUUD state on the frustrated 
square latticci^ However, in contrast to the frustrated 
square lattice, the M/Mgat = 1/3 ground state of the 
SSL has no local continuous degeneracies. Therefore, the 
selection mechanism of the UUD state in the SSL is more 
similar to the triangular lattice i^^'^^ It may be interest- 
ing to note that also the magnetization curve of the Ising 
model on the Shastry-Sutherland lattice exhibits exactly 
one plateau with M/M^at = 1/3,^^ 
Furthermore, we performed Monte-Carlo simulations and 
obtained a phase diagram in the (h, T) plane for the 
particular magnetic coupling ratio J'/J = 1/2. We 
confirmed that the Shastry-Sutherland lattice behaves 
like the triangular lattice and it presents two quasi-long- 
range-ordered phase below and above the long-ranged- 
ordered UUD phase. We also showed that the pseudo- 
plateau survives for a small variation e of the magnetic 
coupling ratio around J'/J = 1/2. The phase dia- 
gram was found to present a new incommensurate um- 
brella phase in the low-temperaturc region, at least for 
J'/J = 0.4. 

In the future, it would be interesting to compute order 
parameters^ and characterize the phase transitions of 
the SSL with classical Heisenberg spins. Furthermore, 
magnetic anisotropies are also expected to be important 
for the rare-earth tetraboridesj ^^i^''i^^'^^i^°'^^i^^ Finally, it 
has been observed on some other lattice o^^i^°'^^ that cou- 
pling to phonons leads to a stabilization of magnetization 
plateaux. It would therefore be interesting to study the 
classical Heisenberg model coupled to lattice degrees of 
freedom also on the Shastry-Sutherland lattice. 
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